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ABSTRACT 



Simultaneous measurements of water surface fluctuations 
and horizontal water particle velocities in a line perpen- 
dicular to the direction of wave approach extending across 
the surf zone were taken in varying surf conditions at 
two locations. The spectral velocities calculated using 
linear theory as a transfer function underestimated measured 
values by 79-86% at the peak of the spectrum. The coherence 
values were generally low indicating non-linear and turbulent 
conditions. Strong harmonics in the spectra of the waves 
and water particle velocities further suggest a non-linear 
system. The theoretical phases computed using linear 
theory did not accurately predict the observed phases. In 
general breaking waves can be characterized as a strongly 
non-linear wave phenomenon. Measured frequency distributions 
were compared with both Gaussian and Gram-Charlier distribu- 
tions by using the chi-square goodness-of-fit test. Quali- 
tatively, the Gram-Charlier distribution gave the better fit 
to the flow velocity data. 
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I. INTRODUCTION 



Wave theories developed for deep water waves can be 
applied with some degree of certainty, and can be tested 
in laboratory and field situations. These theories general- 
ly hold until such time that the shoaling process begins, 
and with somewhat lesser accuracy, throughout the shoaling 
process up to the point of near breaking. At the breaker 
point, however, there is a transition from ordered to 
turbulent motion and the description of wave kinematics 
becomes more difficult. The most forthright approach to 
the problem of describing the kinematics of surf zone 
breaking waves is through direct measurements . Advances in 
instrument design have led to simple, sturdy, measuring 
devices with rapid response time which can sense the small 
scale as well as the large scale features. 

The' study of the kinematics of breaking waves in the 
surf zone has progressed slowly due both to the problems 
encountered in making direct field measurements and the 
difficulty in modeling the surf zone in the laboratory. 
Breaking waves were measured in the laboratory in early 
work by Iverson (1953) using photographic techniques. 

The length of the channel limited the wave type to plunging 
and surging breakers. Adeyemo (1970) made similar labora- 
tory measurements using hydrogen bubble and photographic 
methods. Gaughan and Komar (1975) applied the theory of 
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wave propagation in water of gradually varying depth, as 
developed by Biesel, to determine the dependence of breaker 
type on the beach slope tangent and the deep water wave 
s teepness . 

Inman and Nasu (1956) made field measurements of water 
particle velocity by measuring the drag force under the 
wave in order to infer particle velocities. Miller and 
Ziegler (1964) used both acoustic and electromagnetic 
current meters to determine the particle motion in the 
surf zone. Walker (1969) made studies using propeller- type 
flow meters. Wood (1973) measured waves and currents in 
the surf zone using movies of dye movement and capacitance 
wave gauges. Fiihrboter and Biisching (1974) utilized a two- 
component current meter and two pressure-type wave meters 
to measure simultaneous orbital velocities and water levels. 
Thornton (1968) , Steer (1972) , Thornton and Richardson 
(1973) and Bub (1974) used pressure meters, capacitance 
wave gauges and electromagnetic current meters to measure 
surface profiles and particle velocities; the work 
presented here is an extension of these studies. 

The objective of this research is to study the kine- 
matics of water particles in breaking waves within the 
surf zone. Simultaneous measurements were made of the 
instantaneous wave profile and the horizontal water 
particle velocities at fixed locations in the surf zone. 
Estimates of the probability density functions and powci' 
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spectra of the wave heights and particle velocities are 
made. The applicability of using linear wave theory as 
a spectral transfer function in computing velocity spectral 
components from the power spectrum of the waves is 
measured. The computed velocity spectra are compared with 
actually measured velocity spectra. Theoretical phase 
spectra between wave gauges are calculated and compared 
to measured phase spectra. 
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II. MEASUREMENTS 



A. EXPERIMENT SITES 

The experimental sites are in the vicinity of Monterey, 
California. The beaches here are some of the first inten- 
sively studied to gain an understanding of amphibious war- 
fare techniques and are described by Bascom (1964) . It was 
desired to measure the various types of breaking waves in 
the study including spilling, plunging and surging breakers. 
The manner in which waves break depends very much on the 
characteristics of the beach and near-shore bottom slope 
(Table I). Plunging and spilling breakers occur most 
frequently at the Del Monte Beach site within Monterey Bay, 
The waves at this location are generally severely direc- 
tionally filtered and refracted by the geometry of the bay 
and impinge almost perpendicular to the shore with a 
resulting simplification in the wave description. A second 
experimental site was Carmel River Beach, five miles to 
the south, where the beach is very steep and very often 
has surging type breakers. Again here, the beach is within 
an embayment and the waves impinge almost perpendicular to 
the shore. A typical beach profile and instrument loca- 
tion for Del Monte Beach is shown in Figure 1. Del Monte 
Beach profiles for 4-10 March 1975 are given in 
Appendix A. Figure 2 shows the beach profile and instrument 
location at Carmel River Beach on 29 May 1975. A beach 
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TABLE I. BEACH AND WAVE CHARACTERISTICS 
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profile taken the following day was identical to 
Figure 2. 

B. INSTRUMENTATION 

Many problems were encountered in the initial attempt 
in February 1975 to measure the desired parameters. The 
first wave gauges were manufactured from 1.0 cm diameter 
stripped RG-11 coaxial cable as was done by Bub (1974) . 
Large masses of kelp torn loose by high waves, flotsom 
and the pounding surf easily broke these cables. Scouring 
around the base-plates of the instrument towers and the 
blade-type anchors caused the towers to topple. The con- 
struction of the wave gauges was modified to a three- 
eighths inch outside diameter stainless steel rod, covered 
with a one-sixteenth inch wall thickness teflon tubing 
fitted tightly over the rod; these wave gauges withstood 
the severe forces in the surf well. Steel pipe extensions 
each 0.6 m in length were added under the tower base plates 
and resulted in stable platforms. An extension was built 
for the top of the anchors; with this device, the anchors 
were sledge-hammered deeper into the beach so that the top 
of each anchor was approximately 15 cm below the sand. The 
extension was then removed and scouring was minimized. 

Measurements at Del Monte Beach were made using two 
electromagnetic flow meters and three capacitance wave 
gauges. One flow meter and two wave gauges were used at 
Carmel River Beach. The instruments were mounted on the 



17 



towers within the surf zone. All equipment was calibrated 
in the laboratory prior to the experiments. 

The electronics package for each of the 2.5m 
capacitance wave gauges was constructed using the design 
of McGoldrick (1969) . The capacitance wave gauge operates 
on the principle that a change of capacitor plate dimensions 
changes the circuitry voltage output. The rod and the 
sea water act as plates and the teflon insulation as the 
dielectric. As the wave elevation fluctuates, the capaci- 
tance of the circuit changes and the voltage output responds 
linearly. These water surface fluctuations were sensed by 
a transistorized circuit operated by self-contained bat- 
teries. The electronics package \\ras housed in a water- 
tight brass case which was mounted on the tower. This 
enabled the connecting leads to be less than 30 cm, thereby 
minimizing wire-to-wire capacitance. Accuracy was estimated 
to be ±0.005 m. The calibration plots are shown in 
Appendix B. Minor modifications to the wave gauges 
immediately following the Del Monte Beach measurements 
are the cause for the different calibration graphs in 
Appendix B. 

The two flow meters were Marsh-McBirney Model 721 
Electromagnetic Current Meters. The flow meter operation 
is based on Faraday’s principle of electromagnetic induction. 
Each probe measures water velocity in two orthogonal 
directions. The flow meters were calibrated with an 
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oscillating platform attached to an eccentric arm driven 
by a variable speed motor using the method of Thornton 
and Krapohl (1974) . Measurement accuracy was determined 
to be ±0.005 m/sec during calibration. The flow meter 
calibration plots are shown in Appendix C. 

In the field the instruments were attached to 6.3 cm 
outside diameter steel pipes which were 3.6 m high with a 
1.0m base plate and the 0.6 m steel pipe extensions. The 
towers were placed on a line perpendicular to the predomi- 
nant wave direction and were erected during low tide when 
the beach was easily accessible. The measurements were 
then conducted at high tide. The tidal range in the 
Monterey area is typically two meters. The towers were 
supported by steel guy wires fastened to the anchors. 
Several types of anchors were tried and the blade anchor 
worked quite satisfactorily in the relatively coarse sand. 
The wave gauges were attached to the towers with locally 
fabricated micarta clamps. A 1.6 cm steel stiffening rod 
was added to the arrangement in order to reduce wave gauge 
vibration following breaker impact. The flow meters were 
positioned directly under the wave gauges with the axes 
aligned horizontally ( + X, shoreward) and vertically (+Y, 
upward). A carpenter’s level was used to establish axis 
alignment with an estimated error of ±2 degrees. A typical 
tower and sensor arrangement is shown in Figure 3. 

All signals were cabled ashore and recorded on a Vidar 
Corporation 32-channel digital data acquisition system. 
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FIGURE 3. Schematic of Typical Instrument Tower with Sensors 



A Sangamo model 3500 14 -channel FM tape recorder was 



utilized as 
strip chart 
performance 
appropriate 



a secondary recording system. A Brush 8-channel 
recorder was used to monitor the instrumentation 
during recording and as a means to select the 
data sections to be analyzed. 
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III. ANALYSIS OF DATA 



Record lengths of 30 minutes from each data set \srere 
analyzed. The details of the analysis system including 
calibration factors are given in Appendices D and E. A 
mean value was calculated for all data sets and the data 
was linearly detrended to exclude the rise and fall of the 
tide. The variance, standard deviation and average period 
were calculated. The average period was determined by 
calculating the time between zero upcrossings. However, 
perturbations caused by secondary gravity waves "riding" 
the primary gravity waves, capillary waves and minor instru- 
ment noises increased zero upcrossing occurrences and 
caused the calculated average periods to be lower than 
visually observed. Probability density functions for each 
data set were calculated and graphically compared with 
Gaussian and Gram-Charlier distributions. 

For each signal an auto-covariance function was calcu- 
lated and smoothed with a Parzen window. A Fourier trans- 
form was then applied to the smoothed auto-covariance 
function and the power spectrum determined. A cross- 
covariance function between data sets was computed, smoothed 
with a Parzen window, a Fourier transform applied, and the 
cross - spectrum computed. The coherence and phase were then 
calculated. To compare waves and flow velocity, the wave 
profile spectrum was converted to a theoretical velocity 
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spectrum for comparison with the measured flow velocity 
spectrum. The two power spectra, the coherence and the 
phase were then plotted. In comparing two wave profiles, 
theoretical celerities were calculated using linear theory 
and its shallow water approximation which were then con- 
verted to phases and plotted for comparison with the 
measured phase difference. 

Aliasing testing was performed on all spectra. 

Nyquist frequencies of 1.25 and 0.98 Hz were determined 
optimal for Del Monte Beach and Carmel River Beach data 
sets, respectively. The frequency bands of 0.0 to 1.25 
and 0.98 Hz include the region of gravity waves (0.033 
to 1.0 Hz) . 

The maximum lag time in calculating the covariance 
functions was taken as five percent of the record giving 
a spectral bandwidth resolution of 0.00556 Hz and results 
in 40 degrees of freedom for each spectral estimate. The 
80 percent confidence limits for 40 degrees of freedom 
using a chi-square distribution are found to be between 
0.73 and 1.30 of the measured power spectral estimates. 
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IV. RESULTS 



A. QUALITATIVE DESCRIPTION 

A number of similarities of wave form can be observed 
for various types of breakers occurring on different 
beaches around the world. Figures 4 and 5 are typical 
analog records of waves and horizontal velocities beneath 
the waves obtained from Del Monte and Carmel River Beaches, 
respectively. In general, there is a quick drawdown of 
water just before the breaking wave arrives, followed by 
a steep, almost vertical leading edge, and a sloping pro- 
file toward the trailing edge. On the trailing edge of 
the wave secondary waves are often noted. These are 
harmonics of the primary wave frequency and are indicative 
of very non-linear waves. At Del Monte Beach the waves 
break as plunging and spilling breakers giving a generally 
saw-toothed shape; the plunging or spilling occurs rapidly 
at the crest and moves down the wave. The surging breakers 
occurring at Carmel River Beach rise and fall off more 
gradually . 

B. PROBABILITY DENSITY FUNCTIONS 

1 . Gaussian and Gram-Charlier Frequency Distributions 
Frequency distributions were computed for each 
measured quantity for all runs and compared with Gaussian 
and Gram-Charlier distributions. The chi-square 
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FIGURE 4. Typical Analog Record of Waves and Horizontal Velocities 

beneath the Waves from Del Monte Beach 



I 
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FIGURE 5. Typical Analog Record of Waves and Horizontal Velocities 

beneath the Waves from Carmel River Beach 



goodness -of- fit test was computed using both theoretical 
distributions. The variance, standard deviation, skewness 
and kurtosis were also computed for all data sets and are 
summarized in Table II. The closer the chi-square fit 
parameter is to zero, the better the fit. 

In deep water, the sea surface has approximately a 
Gaussian distribution which results in a skewness of zero 
and a kurtosis of 3. However, non-linearities introduce 
skewness and kurtosis values that deviate from Gaussian 
and result in a distribution more closely approximated by 
the Gram-Charlier . These non-linearities are evident in 
the frequency distribution of Wave Gauge #2 on 4 March as 
shown in Figure 6, and result in bi-modal or tri-modal 
distributions. Figure 7 is the graph of distributions 
computed for the particle velocity at Flow Meter #2 on the 
same date. The latter has a uni-modal distribution. 

Appendix F contains distributions computed for all other 
data sets. 

The distributions for all flow velocities are either 
uni-modal or have less pronounced secondary peaks compared 
to the distributions of wave heights. As noted in Table II, 
this results in nearly equal fit parameters for both Gaussian 
and Gram-Charlier distributions, although the latter values 
are smaller, and therefore give the better fit in ten of 
eleven cases. The exception is Flow Meter #2 on 8 March. 
Additionally, the fit parameters for velocity components 
are much closer to zero than are the values computed for 
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TABLE II. COMPUTED STATISTICS 





d 




















































0 


















































U 


H 


















































0 


! ( 


















































4-> 


d 


CO 


CO 


Cst 


LO 


rH 


n oo 


LO 


LO 


vO 


NT LO 


O 


IN 


or 


IN 


to 


Gh 


Gh 


IN Gh 


LO 


O 


rH 


0 


d 




to 


CO 


(Nt 


CO 


LO 


nT 


o 


rH 


in 


CNI 


oo 


to 


to 


CNI OO 


Oh 


O 


nT 


CNI 


nT 


to 


N 


Gh 


6 


-d 


















































d 


CJ 


iH 


LO 


oo 


OO 


to 


nT 


nT 


NT 


LO 


VO 


in 


OO 


CNI 


CNI 


CNI 


rH 


vO 


CNI 


CNI 


NT 


Gh 


Ni- 


O 


vO 


U 


1 


T— { CO 


IN 


CO 


LO 


vO 




LO 


CO 


vO 


rH 


to 


NT 


O 


O 


Oh 


o 


oo 


vO 


rH 


oo 


ce 


Gh 


nT 


d 


£ 


NT 




oo 




LO 




VO 




NT 




NT 




to 


IN 


rH 


rH 


LO 


to 




VO 




rH 




CNI 




cd 










vO 




rH 




rH 










00 
















rH 








d 










rH 








LO 










rH 






















O 


CO 


















































U 




















































d 




















































3 




















































c r 


d 


















































in 


cd 


<NJ 


to 


CNI 


o 


LO 


CNI 


O 


CO 


to 


CO 


o 


NT 


r^ 


O 


NT 


O 


to 


CNI O 


vO vO 


O 


CO 


CNI 


i 


H 


CM 


LO 


in 


"=3* 


to 


to 


to 


rH 


N 


O CO 


vO 


to 


O 


(NI 


rH 


rH 


Gh 


OO 


IN 


LO 


LO 


NT 


LO 


H 


CO 


















































rj 


CO 


CO 


oo 


CO 


CO 


LO 


CNI 


vO 


LO 


LO 


nT 




O 


in 


CNI 


OO 


LO 


VO 


OO 


to 


O 


00 


oo 


vO 


vO 


O 


d 


CO 


o 


VO 


o 


to 


CO 


in 


oo 


o 


in 


n 


NT 


o 


rH 


to 


O 


O 


LO 


oo 


CNI 


rH 


CO 


Gh 


CO 




cd 


CO 


CNI 


oo 


rH 


oo 




Nt 




o 




! 1 




to 


to 


rH 


nT 


to 


to 




CNI 


rH 


nT 




to 




O 






rH 




CNI 








rH 






















rH 




rH 






CO 


CO 


to 


O 


IN 


to 


o 


CO 


in 


to 


oo 


00 


O 


LO 


to 


O 


vO 


Oh 


o 


oo 


v0 


vO 


rH 


NT 


CNI 


•r 




VO 


to 


oo 


00 


LO 


CO 


NT 


NT 


o 


CNI 


o 


rH 


rH 


IN 


CNI 


nT 


rH 


rH 


vO 


O 


CNI 


rH 


to 


Gh 


CO 


CO 


to 


Cst 


o 


LO 


to 


o 


CO 


vO 


rH 


vO 


oo 


O 


Gh 


to 


to 


O 


NT 


to 


LO 


to 


NT 


rH 


OO 


C 


> 


CO 


rH 


(NJ 


vO 


rH 


Nt" 


LO 


CO 


CO 


CNI 


NT 


NT 


rH 


to 


vO 


o 


NT 


vO 


oo 


to 


LO 


IN 


Gh 


to 


4- 


j 


N 


rH 


to 


LO 


LO 


CO 


in 


CO 


rH 


o 


CO 


o 


to 


to 


LO 


to 


NT 


rH 


nT 


rH 


to 


vO 


LO 


IN 






















































5 


<NJ 


to 


to 


to 


LO 


CNI VO 


CNI 


Nt* 


to 


nT 


to 


to 


oo 


CNI tO 


NT 


NT 


to 


to 


to 


LO 


NT 


to 












to 








rH 
































CO 


vD 


rH 


CO 


LO 


LO 


(Nt vO 


to 


vO 


CO 


in 


to 


nT 


to 


oo 


LO 


Oh 


CNI OO 


to 


(NI 


oo 


vO 


vO 


CO 




rH 


LO 


vD 


to 


oo 


o 


00 


o 


oo 


oo 


CNI 


to 


Gh 


NT Oh 


to 


nT 


to 


CNI 


nT 


N 


OO 


to 


0 


tO 


in 


LO 


to 


nT 


LO 


to 


LO 


00 


CO 


CO 


CNI 


IN 


LO 


N 


IN 


Oh 


to 


N 


Gh 


vO 


CNI 


Gh 


nT 


d 


VO tO 


CO 


IN 


LO 




o 


CO 


CO 


o n 


IN 


rH 


N 


oo 


Oh 


to 


IN 


IN 


vO 


IN 


LO 


to 


NT 






vO 


to 


oo 


(Nt 


CO 


iH 


oo 


rH 


nT 


rH vO 


o 


O 


Gh 


o 


LO 


o 


NT O 


IN 


rH 


LO 


CNI 


vO 


0 


















































HP 


o 


o 


o 


O 


to 


O 


o 


O 


rH 


O 


o 


o 


O 


o 


o 


o 


o 


o 


o 


O 


O 


o 


O 


O 


CO 




1 












1 








I 


1 




















i 






d 


















































"d 


o 


NT 


rH 


LO 


rH 


to 


NT 


o 


CO 


oo 


N 


CO 


LO 


nT 


to 


vO 


Oh 


oo 


nT 


vO 


oo 


vO 


o 


LO 


CNI 


?H ■ 


*H 


N 


iH 


iH 


O 


rH 


rH 


o 


NT 


nT 


LO 


CO 


oo 


vO 


Gh 


IN 


rH 


(NI 


nT 


rH 


Gh vO 


o 


OO 


rH 


d 


4~> 


CO 


to 


o 


(Nt 


CO 


LO 


LO 


to 


o 


rH 


vO 


n 


Gh 


rH 


to 


NT 


IN 


rH 


IN 


NT vO 


Ch 


to 


vO 


“d 


cd 


rH 


VO 


vO 


to 


NT 


o 


vO 


vO 




VO 


CNI 


(NI 


CNI 


(NI 


(NI 


to 


oo 


NT 


CNI 


rH 


Gh 


Gh 


Gh 


Gh 


d 


•H 


(NI 


CO 


rH 


00 


rH 


vO 


rH 


vO 


o 


LO 


rH 


nT 


rH 


rH 


nT 


rH 


LO 


rH 


LO 


rH 


LO 


rH 


vO 


O 


d 


> 


















































4-> 


0 


o 


o 


O 


o 


O 


O 


o 


o 


o 


o 


O 


o 


O 


O 


o o 


o 


O 


O 


O 


o 


O 


O 


O 


CO 


Q 


















































0 


















































u 


LO 


CO 


LO 


nt 


NT 


o 


to 


CNI 




vO 


O 


LO 


rH 


IN 


oo 


rH 


o 


rH 


Gh 


CNI 


rH 


O 


to 


NT 


d 


oo 


LO 


vO 


CNI 


(NI 


CNI 


CNI 


CNI 


CO 


to 


rH 


o 


oo 


oo 


LO O 


Oh 


O 


OO 


CNI 


o 


VO 


NT 


CNI 


d 


IN 


n. 


LO 


(Nt 


CNI 


vO 


n 


O 


nT 


LO 


VO 


to 


VO 


nT 


Oh 


CO 


nT 


O 


N 


to 


vO 


Gh 


rH 


Gh 


•H 


nT 


(NI 


(Nt 


CO 


(Nt 


vO 


CNI 


NT 


o 


rH 


rH 


oo 


rH 


rH 


N 


rH 


nT 


CNI 


IN 


i — 1 


LO 


to 


oo 


O 


U 


o 


CO 


O 


vO 


O 


to 


O 


NT 


o 


to 


O 


rH 


O 


O 


rH 


O 


to 


O 


CNI 


o 


to 


o 


NT 


O 


d 






















































o 


o 


o 


O 


o 


o 


O 


o 


o 


o 


O 


O 


O 


O 


O 


O 


o 


o 


o o 


o 


o 


O 


O 


4- 




















































d 


















































0 


















































6 


rH 


rH 


Cst 


(Nt 


rH 


rH 


CNI 


CNI 


rH 


rH 


CNI 


CNI 


to 


rH 


rH <NI 


CNI 


rH 


rH 


CNI 


CNI 


rH 


rH 


CNI 


3 

f , 


=tt= 


=#= 


=tt= 


=**= 




=#= 


=n= 


=tte 




=t*r 




=#= 


=fc= 


=*te 


=tt= 






=tt= 




=tt= 


=«= 


=tfc 


=tt= 


=tt= 


H 


CO 




CO 


s 


CO 




CO 




CO 




CO 


1 


CO 


CO 




C D 




C0 




CG 




CG 




o 


CO 

r-J 




Ph 




PH 




PH 




PH 




PH ^ Ph 


<• 




PH 




PH 




PH 




PH 




PH 


j2: 


>— i 
1— 1 






















































Pd 
















Pd 










Pd 








Pd 


















u 








u 








o 










u 








u 
















d 








d 








d 










d 








u 








d 






0 


d 








cd 








jd 










cd 








d 














4-* 










I 








































d 












































Gh 






Q 


nt 








LO 








VO 










oo 








Gh 








CNI 







28 



Variance for the wave heights is given in meters squared; variance for the horizontal 
velocity components is given in meters squared per second squared. 

Standard deviation for the wave heights is given in meters; standard deviation for the 
horizontal velocity components is given in meters per second. 
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wave components. This is to be expected since the flow 
meters do not experience the surface irregularities to 
the same degree as do the wave gauges. In seven of 
thirteen fit parameter calculations for wave components, 
the Gaussian distribution provided the better fit. This 
is a result of using a truncated form of the Gram-Charlier 
in calciilating the probability density function, , 
given by 

m_ m. m ^ 

p g C (o = p G cmi.o + <r-H 3 * 2 T» 4 + 7 |_h 6 + ] ci) 

where 

P = Gaussian probability density function, 

VJ 

m^ = skewness of the data record, 
m^ = the kurtosis minus 3, 

H n = Hermite polynomials of degree n. 

Higher order moments could be included in equation (1) 
to possibly improve the fit, but would necessitate analyzing 
longer record lengths in order to maintain confidence. 
Stationary problems would be encountered in the longer 
records, thus producing other errors. Hence, the poor 
fit is a deficiency of the Gram-Charlier distribution in 
surf zone applications. 

2 . Skewness and Kurtosis 

The physical geometry of the waves, which 
characteristically have narrow, steep crests and wide, 
shallow troughs in the surf zone produce a positive 
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skewness, indicating a greater amount of time below, 
rather than above, the mean water level. This is true 
except at Wave Gauge //3 on 6 March. (See Table II and 
Appendix F.) Data from this wave gauge was analyzed on 
this date only, because visual observations noted that 
most waves were actually breaking on this gauge, rather 
than on Wave Gauges #1 and #2 as was otherwise the case. 

It should follow that the wave induced particle 
velocities would have probability density functions similar 
to the waves. However, of the eleven horizontal velocity 
data records, four did not have positive skewness. These 
deviations occurred at Flow Meter #1 on 4 March and 
29 May and Flow Meter # 2 on 5 and 6 March. Reflected 
waves were originally thought to be a feasible explanation 
for the negative skewness. If this were so, Flow Meter #2 
on 4 March and 29 May would also expect to have negative 
values. This was not the case. The reason for the large 
skewness and kurtosis values for Wave Gauge #1 on 5, 6 and 
8 March is not known . Similar anomalous results have been 
found in deeper water outside the surf zone by Thornton 
and Krapohl (1974) . 

Kurtosis indicates the peakedness of each 
parameter. Just prior to breaking the waves achieve the 
greatest degree of peakedness and should have the highest 
kurtosis value at this time. Visual observations noted 
that most waves were breaking at Wave Gauge //I on 
4 March and 29 May, at Wave Gauge 112 on 5 March and on 
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Wave Gauge #3 on 6 March. It should be expected that the 
kurtosis values would be greatest at these given gauges 
for the respective dates. With the exception of 29 May, 
this did not turn out to be true. The wide surf zone 
caused by varying breakers heights on 8 and 9 March make 
evaluation of kurtosis values difficult on these dates. 

C. COMPARISONS OF THEORETICAL AND MEASURED POWER SPECTRA 
The theoretical velocity spectra generated from the 
wave spectra and the measured velocity spectra were 
examined in order to determine the applicability of linear 
(Airy) wave theory in the breaker zone. Thornton and 
Richardson (1974) measured waves with maximum heights of 
1 m and showed that the application of a linear theory 
transfer function to the surface profile spectrum resulted 
in a calculation of the horizontal water particle velocity 
that underpredicted measured wave-induced horizontal 
velocity spectral components by about 50 percent. The 
coherence values between waves and horizontal velocity 
were high, ranging above 0.75 for an approximate frequency 
range of 0.075 to 0.60 Hz. The phase angle computations 
showed the calculated velocity components leading the 
measured velocity components by an average of 20 degrees. 
Bub (1974) took measurements in mild surf resulting in 
theoretical horizontal velocity values 13 percent lower 
than measured velocity. Coherence values varied from 
0.5 to nearly 1.0 over a frequency range of 0.10 to 0.65 Hz 
An approximate zero degree phase was observed out to 1.0 Hz 
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For this study wave heights of approximately 1 m 
were typically observed at Del Monte Beach and up to 2 m 
were recorded at Carmel River Beach. The theoretical 
horizontal velocity spectra were computed from the wave 
spectra using the linear theory transfer function 



uc 



w - t p v f > 



( 2 ) 



where 

S uc (f) = calculated horizontal velocity power 

? 

spectrum (m /sec) , 

S n (f) = measured wave height power spectrum (m -sec), 
o = angular frequency = 27rf(sec 1 ) , 
k = wave number (m ^) , 
h = mean water depth (m) , 

z = depth of flow meter below the mean water 
depth (m) . 

Figures 8 and 9 are graphs of the measured and calculated 
spectral velocity values at Tower #1 for 4 March and 
29 May. Additional spectra are contained in Appendix G. 

On the average at Del Monte Beach, the calculated hori- 
zontal velocity values underestimated the measured veloci- 
ties by 79%. The same theoretical calculations for the 
data collected at Carmel River Beach underestimated the 
measured values by 86%. 

A possible explanation for the low calculated values 
is that the wave profiles change during the shoaling and 
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FIGURE 8. Power, Coherence and Phase Spectra for Flow 
Meter HI and Wave Gauge If 1 on 4 March 197 5 
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Meter #1 and Wave Gauge Hi on 2 9 May 1975 
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breaking process, becoming more peaked at the crest, 
thus lessening the applicability of linear (Airy) theory 
which assumes sinusoidal wave shapes. Additionally, the 
wider separation between the measured and theoretical 
spectral values above 0.5 Hz may be explained by turbu- 
lence of the breaking waves in the higher frequency areas. 

D. EXAMINATION OF SPECTRAL HARMONICS 

An examination of the low frequency spectral components 
for each data set indicates that pronounced approximate 
harmonics are present in the data analyzed from Del Monte 
Beach. Such harmonics are not evident in the Carmel River 
Beach spectra. The frequencies and amplitudes of each 
analyzed spectral peak is listed in Table III. The fre- 
quency differences between spectral peaks and the ratios 
of the second, third and fourth peaks versus the primary 
peak are listed in Table IV. 

Spectral analysis shows that definite second, third 
and fourth harmonics are present for 4 March and 5 March 
with few exceptions. Evidence of second and third har- 
monics is present in the analyzed data from Wage Gauges 
#1 and #2 on 6, 8 and 9 March. An examination of the 
spectral amplitudes shows that there is no consistent 
relationship between various components either for a 
given data set or when comparing like components of 
different data sets. 

The maximum wave energies were concentrated at frequen- 
cies ranging from 0.061 Hz to 0.092 1 1 z corresponding to 
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TABLE III. FREQUENCIES AND AMPLITUDES AT SPECTRAL PEAKS 
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£ = frequency (Hz) at spectral peaks 

a = amplitude at spectral peaks 

subscripts: 0 = primary, 1 = secondary, etc. 



TABLE IV. HARMONIC FREQUENCIES AND SPECTRAL AMPLITUDE RATIOS 



o 


vO 


[>. 


00 


rH 


CP 


oo 


OO 


CM 




























H 


LO 




05 


CM 


rH 


CM 


to 




CM 


LO 


to 


i 


i 


i 


1 


I 


I 


1 


i 


1 


i 


1 


1 


i 


OO 


CP 


1 


\ 


o 


O 


o 


o 


O 


o 


o 


o 


I 


i 


i 


1 


l 


l 


1 


i 


1 


i 


1 


1 


i 


H 


o 


i 


to 












































• 


• 






















































o 


CM 


CP 


o 


VO 


o 


CM 


CM 


o 


LO 


o 




o 


CP 


OO 


cp 


CM 


H" 


CP 


CM 


00 


to 




to 


o 


oj 


LO 


to 


LO 


LO 


H 


LO 


CP 


to 




00 H- 


LO 


to 


00 


to 


LO 


LO 


H 


CP 


IN- 


to 


OO 


CM 


CM 




O 


o 


o 


o 


H 


o 


O 


o 


o 


O 


O 


o 


o 


o 


o 


o 


o 


H 


o 


O 


H 


CM 


H 


H 


CM 


















































o5 


















































o 




CO 


cp 


CP 


LO 


oo 


CM VO 


o 




H 


CM 


vO 


vO 


CM 


o 


oo 


O 


CM 


CM 


CM 


o o 


O 


o5 


VO 


CO 


to 


LO 


H 


LO 


CP 


H 


LO 


cp 


H 


o 


LO 


c^ 


o 


CM 


CM 


LO 


oo 


CM 


CP 


CM 


H 


o 




rH 


rH 


rH 


rH 


to 


H 


H 


H 


CM 


H 


H 


H 


o 


H 


H 


H 


H 


CM 


H 


CM 


CM 


H 


LO 


CM 


H 


















































05 












































H 






CM 


rH 


rH 


CM 


rH 




CM 


H 


LO 




























CP 


H 




■U 


VO 


VO 


VO 


LO 


vO vO vO 


H" 


l 


l 


1 


l 


1 


1 


1 


1 


1 


1 


I 


1 


l 


OO 


LO 


l 


l 


O 


O 


O 


o 


O 


O 


O 


O 


I 


i 


I 


l 


1 


1 


1 


1 


1 


t 


I 


1 


1 


o 


O 


l 


tO 












































• 


• 




4-1 


















































rH 


Cxi 


CO 


rH 


CM 


VO 


P"- VO 


CO 


00 


o OO 


H 


N- 


to 


tO H" 


CM 


LO LO 


CP 


rH 


CP 


o 


VO 


4-1 


VO 


vO vO 


C- 


LO 


VO 


LO 




c^ 


o 


N- 


O VO 


00 


OO 


OO 


CP 


oo 


00 


o 


O 


r- 


o 


LO 


i 


o 


o 


O 


o 


o 


o 


o 


o 


o 


H 


O 


rH 


o 


o 


o 


o 


o 


o 


o 


H 


H 


o 


rH 


H 


CM 


















































4-1 


















































O 


rH 


C-. 


CM 


CM 


c- 


H 




H 


oo 


to 


OO 


CO 


cp 


H" 


H" LO 


^3" 


o 


o 


CP 


LO 


H- O 


LO 


4H 


vO VO 


vO 


vO 


VO 


VO 


vO vO 


c^ 








oo 


oo 


CO 


IN. 


00 


o 


o 


O 




*3- 


LO 


CP 


i 


O 


o 


o 


O 


o 


o 


o 


O 


o 


o 


o 


O 


o 


o 


o 


o 


o 


H 


H 


\ — 1 


o 


O 


o 


o 


r—P 


















































4^ 






















































rH 


rH 


rH 


Cx. 


l>~ 




C^ 


oo 


oo 


oo 


CO 


oo 






LO 


CM 


CM 


CM 


CM 


oo 


00 


oo 


o 


vO 


vO vO vO 


VO 


VO 


VO 


VO 


c- 




c- 




c- 


00 


oo oo 




CP 


CP 


CP 


CP 


IN. 


N. 


CM 


4H 


O 


o 


O 


O 


o 


o 


o 


O 


o 


o 


o 


o 


o 


o 


o 


o o 


o o 


o 


o 


o 


o 


H 




















































C 


















































0 


















































e 


rH 


rH 


CM 


CM 


H 


H 


CM 


CM 


H 


H 


CM 


CM 


to 


H 


H 


CM 


CM 


H 


H 


CM 


CM 


H 


H 


CM 


3 




=tt= 


=tt= 


=**= 


=tt= 






=tfc 


=tfc 


=te 


=tt= 


=tt= 


=4= 




=3= 


=4t= 


=*fc 


=tt= 


=4= 


=tte 


=tt= 


=tt= 




=*t= 


H 
























































CO 


s 


CO 




CO 








CO 


CO 


S 


CO 




CO 


s 


co 




CO 




CO 


(/) 




IX, 




Ph 




PP 


V?- 


PP 




PP 




PP 






pp 




pp 




pp 




PP 




PP 




C 


















































— t 




















































rC 
















x: 










rC 








& 
















U 








U 








u 










CJ 








o 








>X 








H 








U 








u 










P< 






















0 










05 








o5 










05 








iS 








Sp 








S 
















































05 












































CP 






Q 


H- 








LO 








VO 










oo 








CP 








CM 







39 



£ = frequency (Hz) at spectral peaks 

a = amplitude at spectral peaks 

Subscripts: 0 = primary, 1 = secondary, etc. 



periods of 16.4 to 10.9 seconds, respectively, which agreed 
with visual observations. Additionally, prominent sub- 
harmonic peaks between 0.011 and 0.022 Hz are noted in the 
measured velocity spectra computed from Del Monte Beach 
data. This agrees with the work of Bub (1974), who 
observed a sub-harmonic at 0.011 Hz. Guzza and Davis (1974) 
state that in a theoretical analysis of edge waves excited 
by incoming waves, the prominent edge wave mode is the 
first sub-harmonic of the primary frequency. It is plausi- 
ble, therefore, that the observed low frequency sub-harmonic 
peaks were caused by edge waves in the surf zone. Such 
sub -harmonics are not observed in the measured velocity 
spectra generated from Carmel River Beach data. 

E. PHASE AND COHERENCE 

The curling crests of the unstable breaking waves 
tend to lead the horizontal particle velocities in the body 
of the waves as can be seen in the phase measurements 
between waves and velocities. (See Table V.) The maximum 
phase difference would be expected at breaking. Data from 
Tower #1 on 4 March indicates that the theoretical hori- 
zontal velocity leads the measured horizontal velocity by 
about 20 degrees. Since most waves observed were breaking 
on Tower #1, a possible explanation for the 30 degree phase 
angle at Tower 2 is that the waves reorganized between the 
toilers and broke a second time at Tower #2. Since this 
author’s observations concentrated primarily on the initial 
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TABLE V. PHASE RELATIONSHIPS 
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29 May WG #1 - FM #1 10° to 0.58 Hz At Tower #1 



breaking position, the possible reorganization and second 
breaking were not observed. The phase angles on 5 March 
are 5 and 20 degrees at Towers #1 and #2, respectively. 

This agrees with observations and indicates that the 
breakers started to curl at Tower #1 and broke at Tower #2. 
The zero degree phase relationships on 6 March show that 
the waves had not started to break until shoreward of 
Tower #2. This agrees with the observation of waves 
breaking on Tower #3. The phase angles on 8 and 9 March 
indicate that the smaller breakers began to curl at Tower #1 
and were near breaking at Tower #2. Again, this agrees 
with observations. 

In contrast to the plunging- type breakers at Del Monte 
Beach, surging- spilling breakers were observed at Carmel 
River Beach. This was to be expected considering the beach 
slopes. A spilling breaker is characterized by turbulent 
water forming at the wave crest and eventually flowing down 
and covering the leading edge of the wave. Spilling begins 
at the crest when a small tongue of water moves faster 
than the wave form as a whole (Galvin, 1972). The measured 
phase angle at Carmel River Beach shows that the cascading 
crests led the horizontal particle velocities in the body 
of the wave by 10 degrees. 

The values of coherence show varying degrees of linearity 
in the transfer functions between different spectra. Co- 
herence between wave and horizontal velocity spectra at 
Tower #1 on 29 May is above 0.53 out to 0.58 Hz (See 



4 2 



Figure 9), above 0.68 out to 0.53 Hz on 8 March and above 
0.60 out to 0.41 Hz on 9 March. (See Appendix G.) Other 
coherence calculations show a fair to poor degree of 
linearity, which is consistent with the observed strong 
non- linearit ies . 

F. THEORETICAL AND CALCULATED PHASE SPECTRA 

A theoretical phase spectrum was calculated for com- 
parison with the measured wave spectrum from two wave 
gauges. The theoretical phase was computed from the 
relationship 

= S 2 - (3) 

L 1 

where 

cf> = phase, 

a = angular frequency = 2ir£ (sec , 
x = distance between wave gauges (m) , 

1/7 

C^ = celerity = (gh) ' (m/sec) , 

2 

g = gravity (m/sec ), 
h = mean water depth (m) . 

Figure 10 shows typical examples of the theoretical and 
measured phases. The theoretical phases computed from Del 
Monte Beach data underpredict the actual phase relationships, 
whereas an overprediction is noted in the Carmel River Beach 
phases. The author then determined what value of celerity 
(C 2 D would cause the theoretical and actual phases to 
agree out to the first zero crossing (where <j> = 180° = ir) 
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PHRSE PHRSE 

-ISO 0 iSO -ISO 0 ISO 




- - - THEORETICAL PHASE 
MEASURED PHASE 




FIGURE 10. Examples of Theoretical and Measured Phases 
from Del Monte Beach (top) and Carmel River 

Beach 
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by first solving for from an iteration of equation (3) 

C 2 = 2fx (4) 

and then solving 

C 2 = (g(h+a)) 1//2 (5) 

for a. The results are shown in Table VI. No evident 
relationship between a and other parameters, such as beach 
slope, mean w r ater depth and standard deviation, could be 
determined. 



TABLE VI. PARAMETERS CALCULATED FOR PHASE COMPARISONS 
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V. CONCLUSIONS 



Breaking waves in the surf zone can be characterized 
as being highly non-linear. Qualitative observations of 
wave and velocity profiles show secondary waves indicative 
of non-linear waves. The probability density functions 
calculated for the flow velocity records compared better 
with the Gram- Chari ier distribution than with the Gaussian 
distribution when tested using the chi-square goodness-of- 
fit test. However, seven of the thirteen wave records 
more closely approximated the Gaussian distribution. This 
is a result of the truncated form of the Gram-Charlier 
probability density function used in calculations, and 
points out the importance of including higher order moments 
in describing wave phenomena in very shallow water. 

The values of the horizontal power spectral components 
calculated from wave spectra using linear theory indicate 
a qualitative, but not a quantitative, relationship. Com- 
paring the results of Bub (1974) , Thornton and Richardson 
(1974) and this author, it is evident that increased 
breaker heights result in larger estimation errors. 
Additional investigation is warranted to determine if 
specific breaker heights and/or types result in specific 
underestimation values. If so, linear theory may be 
applied to achieve a reasonably accurate determination of 
flow velocities. The coherence values between waves and 
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horizontal velocity decrease as breaker height increases, 
indicating that the wave motion becomes more non-linear 
and more turbulent under this condition. 

The theoretical phase calculations based on the linear 
theory approximation for celerity did not accurately 
predict the observed phases in this research. The values 
of a, determined so as to force the theoretical phases to 
agree with the observed values at the first zero crossing, 
show no consistency. 

Use of electromagnetic flow meters and the improved 
capacitance wave gauges and instrument towers permits the 
gathering of accurate continuous data with sturdy and 
reliable equipment. 
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4 MARCH 1975 



APPENDIX A 

BEACH PROFILES AT DEL MONTE BEACH 4-10 MARCH 1975 
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DISTANCE IN METERS 



OUTPUT VOLTAGE (VOLTS) OUTPUT VOLTAGE (VOLTS) 




— f — — — j ” i ' 

O.S 1.0 1.5 2.0 2.5 

DEPTH OF IMMERSION (METERS) 
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3.5-" 

3.0" 

2.5" 

2 . 0 " 



x: wave gauge 1 
□ : wave 9 ou S e ^ 



APPENDIX B 

CAPACITANCE WAVE GAUGE CALIBRATIONS 



x = wave 
d i v^avc 
o : wave 



gauge 1 
gauge 2 
gauge 3 



OUTPUT VOLTAGE (VOLTS) OUTPUT VOLTAGE (VOLTS) 
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APPENDIX C 



FLOW METER CALIBRATIONS 
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APPENDIX D 



ANALYSIS DETAILS AND FLOW CHART 



A complete flow chart of data acquisition and analysis 
is presented on the following page. The parameters used 
during spectral analysis were chosen to balance the longest 
record that the computer could reasonably analyze with the 
best resolution over the frequency range of interest. 
Additionally, consideration was given to the computer run 
time . 

In order to determine the optimum Nyquist frequency for 
data analysis the following equations were used: 

At = (VDT) (NSKIP) (NCHAN) (1) 



and 



f N 2 (At) 



( 2 ) 



where 

At = sampling interval, 

VDT = sampling rate, 

NSKIP = number of samples skipped in 
initial data array, 

NCHAN = number of channels, 
f^ = Nyquist frequency. 

By varying the parameters of equation (1) various Nyquist 
frequencies were obtained and used during preliminary IBM 
360 computer analysis. 
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APPENDIX E 



CALIBRATION FACTORS 



Date 


Ins trument 


GALA* 

(meters) 


CALM** 

(meters/volt) 


4 March 


WG #1 


-1.312 


1.23 


and 


WG #2 


-1.870 


4.19 


5 March 


WG #3 


-1.800 


1.14 


6 March 


WG #1 


-1.422 


1.23 




WG #2 


-1.620 


4.19 




WG #3 


-1.800 


1.14 


8 March 


WG #1 


-1.332 


1.23 




WG #2 


-1.570 


4.19 




WG ft 3 


-1.810 


1.14 


9 March 


WG #1 


-1.492 


1.23 




WG ft 2 


-1.510 


4.19 




WG ft 3 


-1.830 


1.14 


29 May 


WG #1 


-0.028 


0.69 




WG #2 


-0.370 


0.66 


All dates 


FM #1 


0.0 


3.05 




FM ft 2 


0.0 


3.21 


* Calibration 
** Calibration 


Additive Factor 
Multiplicative 


Factor 
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APPENDIX F 



PROBABILITY DENSITY FUNCTIONS VS. 
GRAM-CHARLIER FREQUENCY DISTRIBUTIONS 
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